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ABSTRACT 

Scintillation noise due to the Earth’s turbulent atmosphere can be a dominant noise 
source in high-precision astronomical photometry when observing bright targets from 
the ground. Here we describe the phenomenon of scintillation from its physical origins 
to its effect on photometry. We show that Young’s (1967) scintillation-noise approx¬ 
imation used by many astronomers tends to underestimate the median scintillation 
noise at several major observatories around the world. We show that using median 
atmospheric optical turbulence profiles, which are now available for most sites, pro¬ 
vides a better estimate of the expected scintillation noise and that real-time turbulence 
profiles can be used to precisely characterise the scintillation noise component of con¬ 
temporaneous photometric measurements. This will enable a better understanding 
and calibration of photometric noise sources and the effectiveness of scintillation cor¬ 
rection techniques. We also provide new equations for calculating scintillation noise, 
including for extremely large telescopes where the scintillation noise will actually be 
lower than previously thought. These equations highlight the fact that scintillation 
noise and shot noise have the same dependence on exposure time and so if an observa¬ 
tion is scintillation limited, it will be scintillation limited for all exposure times. The 
ratio of scintillation noise to shot noise is also only weakly dependent on telescope 
diameter and so a bigger telescope may not yield a reduction in fractional scintillation 
noise. 

Key words: planets and satellites: detection - atmospheric effects - instrumentation: 
photometers - methods: observational - site testing - techniques: photometric 


1 INTRODUCTION 


High-precision photometry is key to several branches of as¬ 
tronomical research, including (but not limited to) the study 
of extrasolar planets, astroseismology and the detection of 
small Kuiper-belt objects within our Solar System. The dif¬ 
ficulty with such observations is that, although the targets 
are bright, the variations one needs to detect are often small 
(typically ~0.01% to ~0.1%). This is within the capabili¬ 
ties of modern detectors. However, when the light from the 
star passes through the Earth’s atmosphere, regions of tur¬ 
bulence cause intensity fluctuations (seen as twinkling by 
the naked eye) called scintillation. This scintillation, which 
induces photometric variations in the range of ~0.1% to 
1.0%, limits the detection capabilities of ground based tele¬ 


scopes (e.g. Brovmfc Gillilandl 1 19941 : iHeaslev et al.l 1 19961 : 


I Ryan fc Sandleill 19981).' 


Knowing the level of scintillation noise is important 
because it will enable performance assessment, calibration 
and optimisation of photometric instrumentation. It will 
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also help to explain and constrain model fits to photo¬ 
metric data (for example, e xtrasolar planet transit/eclipse 
light curves; iFohring |2014| ). and to help develop scintil¬ 
lation co rrection concepts su ch as Conjugate-Plane Pho¬ 
tometry dOsborn et al.l l2 01ll ). Tomographic wavefront re¬ 
construction (IOsbornll20i5ir and active deformable mirror 
techniques (IViotto et al.ll201^ . It would also enable passive 
techniques such as ‘lucky photometry’ where only data taken 
during photometric conditions (i.e. in times of low scintilla¬ 
tion noise) are used in the reduction process. 


Young proposed an equation which can be used to es¬ 
timate the scintillation noise for an observation given the 
telescope’s altitude and diameter, and the observation’s ex¬ 
posure time and airmass. This e quation is regularly used by 
many astronomers (for example ISouthworth et al.ll2009l ) to 
estimate the scintillation noise in t heir measurements. Re¬ 
cent work bv iKornilov et ahl (l2012l) showed that this equa¬ 
tion tends to underestimate the median scintillation noise 
by a mean factor of 1.5. 


As well as presenting new results on scintillation, we 
hope this paper will serve as a useful guide for astronomers 
to understand, estimate the size of, and correct for scintilla- 
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tion. We show that using the theoretical scintillation noise 
calculated from the median optical turbulence profile for a 
particular site is a better estimate of the median scintilla¬ 
tion noise. However, we also show that the measured scin¬ 
tillation noise has large variability due to changes in the 
optical turbulence profile. Therefore, if one wishes to know 
the scintillation noise during a particular observation, the 
contemporaneous optical turbulence profile should be used 
if possible. With the proliferation of Adaptive Optics (AO) 
systems, real-time turbulence profilers are becoming more 
common-place at major observatories (see section ITB . As 
the scintillation noise is predominantly caused by high alti¬ 
tude turbulence we can assume that a single profiler at an 
observatory site can be used to precisely predict the scin¬ 
tillation noise applicable to all of the local telescopes. This 
dominance of the high altitude turbulence also explains why 
it is possible to observe highly photometric conditions dur¬ 
ing periods of bad seeing. 

In section [2] the physical origin of scintillation noise 
is described and simulated light curves are shown. Sec¬ 
tion ED reviews Young’s scintillation approximation and in 
section we compare the scintillation noise from Young’s 
approximation with scintillation noise measurements from 
atmospheric profiling instruments. In section 0] we look at 
the magnitude of the scintillation noise on astronomical pho¬ 
tometric measurements. In section [5] we present modified 
equations for the case of extremely large telescopes, when 
the outer scale of the turbulence becomes critical and, in 
section [6l telescopes with large secondary mirrors. 


2 THEORY 

Optical turbulence is caused by the mechanical mixing of 
layers of air with different temperatures and hence density. 
The refractive index of air depends on its density and so tur¬ 
bulence at a boundary of air masses with different tempera¬ 
tures creates a continuous screen of spatially and temporally 
varying refractive indices. 

The wavefront from an astronomical source can be con¬ 
sidered flat at the top of the Earth’s atmosphere. As it prop¬ 
agates to the ground it becomes aberrated by the optical 
turbulence which forms a limit to the precision of measure- 
men ts from ground-bas ed telescopes (for example, see figure 
1 in lOsborn et al.ll201ll l. 

The effect of optical turbulence is twofold. The first ef¬ 
fect is to deform the wavefront by retarding the sections 
passing through regions of higher refractive index. This lim¬ 
its the angular resolution of ground-based telescopes and is 
a first order effect as it depends on the first derivative (i.e. 
local tilt) of the wavefront. The second effect of the turbu¬ 
lence is to locally focus and defocus the wavefront, which 
results in spatial intensity fluctuations, or speckles, in the 
pupil plane of a telescope (figure [T]). This is known as scin¬ 
tillation. Scintillation is a second order effect as it is caused 
by the second derivative or curvature of the wavefront. It 
is the propagation of this curvature that leads to the inten¬ 
sity speckles, therefore it is the higher altitude turbulence 
that is primarily responsible for this effect. This is different 
to the phase aberrations causing images to blur, which is 
dominat ed by the strongest turbulent layer, often near the 
ground (lOsborn et al.ll201Cll l. Therefore, atmospheric seeing 



Figure 1. Example simulated pupil image for a 2.5 m telescope 
and a single turbulent layer at 10 km (left) and an example of a 
real pupil image from the 2.5 m Isaac Newton Telescope on La 
Palma (right). 


(a measure of image quality) and photometric quality can 
be independent. It is possible to observe conditions which 
have bad seeing but good photometric quality. 

The characteristic size of the scintillation speckles is 
given by the radius of the first Fresnel zone, rp = V^iA, 
where 2 is the propagation distance from the turbulent layer 
and A is the wavelength of the light. As a wavefront prop¬ 
agates away from a turbulent layer, increasing 2 and hence 
rp, the spatial intensity fluctuations become larger in terms 
of spatial extent. This is not dependent on the strength of 
the layer, which only affects the magnitude of the intensity 
fluctuations and not their spatial properties; for example, 
a turbulent layer at 10 km observed at 500 nm results in 
speckles of size ~0.07 m, irrespective of the strength of the 
turbulence. 

The intensity speckles traverse the pupil with a velocity 
determined by the wind velocity of the turbulent layer. We 
assume that the speckle evolution timescale is longer than 
the wind crossing time. Speckles from different layers move 
independently and superimpose in the pupil plane. As the 
regions of higher intensity enter and exit the pupil the inte¬ 
grated total intensity also varies. The faster the wind speed 
the quicker the intensity will vary, demonstrating the impor¬ 
tance of the wind speed in estimating the scintillation noise. 
It is these variations which lead to scintillation noise, which 
can limit the precision of photometric measurements. 

The amount of scintillation is referred to as the scin¬ 
tillation index, aj, and is expressed as the variance of the 
normalised intensity from the astronomical target. Assum¬ 
ing no other noise sources, 


2 

' {I? 


( 1 ) 


where / is the measured intensity from the object as a func¬ 
tion of time and () denotes an ensemble average. The scin¬ 
tillation noise is therefore defined as and is the nor¬ 

malised scintillation noise. 

Monte-Carlo computer simulations allow us to study 
the structure and effect of atmospheric scintillation inde¬ 
pendently of any other phenomena. Figure 0] shows a 5 s 
section of a simulated scintillation light curve from an 8 m 
telescope. In this light curve we see noise covering a wide 
range of temporal frequencies, /. Figure |3] shows the power 
spectrum of the scintillation light curve. The power spec¬ 
trum, 'l>x(/)i shows that scintillation noise can be split into 
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Figure 2. Example section of a simulated scintillation light curve. 
The intensity variations here are due entirely to atmospheric scin¬ 
tillation. In this case the telescope has a diameter of 8 m and the 
atmospheric turbulence is given by the mean profile for Paranal. 



Figure 3. Scintillation power spectrum of the light curve shown 
in figure [2] The power spectrum shows significant scintillation 
power over a wide range of frequencies. 


two regimes. Above a certain frequency, the power spec¬ 
trum shows a power-law dependence of -11/3. Below this 
frequency the spectrum is fiat. This turning point is deter¬ 
mined by the amount of spatial averaging of the scintillation 
and is dependent on the telescope diameter and the wind ve¬ 
locity. 


3 ESTIMATING SCINTILLATION NOISE 

3.1 Estimating scintillation noise with Young’s 
approximation 

A standard and often implemented equation for calculating 
the expected scintillation variance at any particular site is 
Yo ung’s approxm ation. This approximation, first suggested 
by lYoung (Il967h . is intended as an approximation and not 


a precise prediction. Young’s approximation is given by 

cry = 10 X 10~®D~'’^®t~^(cos7)“^ exp (—2/iobs/i?), (2) 

where D is the diameter of the telescope, t is the exposure 
time of the observation, 7 is the zenith distance, hobs is the 
altitude of the observatory and H the scale height of the 
atmospheric turbulence, which is generally accepted to be 
approximately 8000 m. All parameters are in standard SI 
units. 

3.2 Estimating scintillation noise with 
Atmospheric Turbulence Profiles 

Young’s equation is an empirical approximation based on 
observations and is not intended to provide a precise predic¬ 
tion. Therefore, a more precise option is to use concurrent 
measurements of the atmospheric turbulence to estimate the 
scintillation noise. 

As the scintillation noise is caused by the intensity 
speckles entering and leaving the telescope pupil there are 
three regimes of scintillation noise: 

• For short exposures the intensity speckles appear frozen 
in the pupil and no temporal averaging occurs. 

• For longer exposures the speckles traverse the pupil dur¬ 
ing an exposure and the scintillation noise will be reduced by 
temporal averaging. The amount by which the scintillation 
noise is reduced is dependent on the wind speed. 

• For small telescopes of approximately the same size as 
the speckles and smaller there will be a significant wave¬ 
length dependence. This is because the size of the speckles, 
rp, IS wavelength dependent. 

The scint i llation index for short exposures is given by 
iKenvon et al.l (l2006l i (see appendix 1X1 for details), 

poo 

erf,SB = 17.34L>"'^/^ (C 0 S 7 )"® / (h) dh, (3) 

Jo 

and for long exposures, 

crj,ie = 10.66D t (cos 7 ) J V±{h) 

where C^{h) is the profile of the refractive index structure 
parameter, a measure of the optical turbulence strength, h 
is the altitude of the turbulent layer, with h = 2 cos ( 7 ), 
V±{h) is the wind velocity profile and a is exponent of the 
airmass and is usually taken to be -3.5. Note that a will 
depend on the wind direction and vary between (cos 7 )”^ 
for the case when the wind is transverse to the azimuthal 
angle of the star, and (cos 7 )”^ in the case of a longitudinal 
wind direction. This difference comes from geometry; if the 
wind direction is parallel to the azimuth of the star, the pro¬ 
jected pupil onto a horizontal layer is stretched by a factor 
of 1 /cos 7 . 

We can re-arrange and solve equations [3] and |4] for ex¬ 
posure time to find an expression for fknee, the exposure 
time at which the scintillation noise moves from the short 
exposure to long exposure regime, 

/■OO I 

iknee = 0.62 D(cOS7)“'^^ / i^) 

Jo y±\'T') 

where —l^a-|-3^0 depending on the wind direction be¬ 
ing parallel or perpendicular to the azimuthal angle of the 
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target, respectively. This change from the short exposure 
regime to the long exposure regime can be seen in the scin¬ 
tillation power spectrum shown in figure [3] For atmospheric 
scintillation the change of regime occurs at an exposure time 
of approximately a few hundredths to a few tenths of a sec¬ 
ond. 

It is interesting to note that in both the short-exposure 
and long-exposure regimes, scintillation noise is independent 
of wavelength. Therefore, in situations where the observa¬ 
tions are scintillation dominated and one is o bserving with 
a mu lti-band instrument, e.g. ULTRACAM dPhillon et al.l 
l 200 'iil . correcting one wavelength channel by another can 
reduce the scintillation noise, assuming that the variability 
that one wishes to measure is restricted to one band. How¬ 
ever, although the scintillation noise is independent of wave¬ 
length, as observations move away from zenith the scintilla¬ 
tion signals become temporally separ ated due to the chro - 
matic dispersion in the atmosphere (IPravins et al. lE^. 
This effect is exacerbated when the wind velocity is aligned 
with the azimuthal angle of the target. In addition, if the 
observation is not scintillation dominated this action can be 
detrimental to the total photometric noise as the other noise 
sources will add in quadrature. 

For small apertures, where the aperture size is smaller 
than the spatial scale of the intensity fluctuations {D < 
rp, i.e. D is smaller than a few tens of centimetres) the 
scintillat ion index can be app roximated using the equation 
given bv IPravins et al.l (I1998 I l 


cTj = 19.2A (cos 7 ) 


70 


Several atmospheric optical turbulence profiling instru¬ 
ments are now in regular use at many of the world’s 
premier observing sites. It is therefore possible to esti¬ 
mate the scintillation noise using tec hniques such as MASS 
(Mul ti Aperture Scintillation System. iTokovinin fc Kornilov! 
2007f), SCIPAR (scin tillation Petection And Ranging, 


Vernin fc R.odc 122^}_and SLOPAR (SLOpe Petection 


And Ranging, WilsonI [2002h . However, only the SLOPAR 
and SCIPAR methods can also measure the profile of the 
wind velocity as a function of altitude, as required by equa¬ 
tion ID 

Here we use stereo-SCIPAR (|SheDherd et al.l 1201,31 : 
lOsborn et al.|[ 201 ^ . a new high-resolution SCIPAR instru¬ 
ment, which can automatically obtain concurrent turbulence 
velocity and turbulence strength profiles. Stereo-SCIPAR 
was installed on the 2.5 m Isaac Newton Telescope (INT), 
La Palma, for three two-week campaigns in 2014 and on the 
1 m Jacobus Kapteyn Telescope (JKT) for three two-week 
campaigns in 2013. An example atmospheric optical turbu¬ 
lence profile for a whole INT night is shown in figure |3] 


3.3 Comparing measured scintillation noise with 
Young’s approximation 

Figure [S] demonstrates the difference between the long ex¬ 
posure scintillation noise as estimated by Young’s approxi¬ 
mation (equation[2j and that calculated from observed SCI¬ 
PAR profiles (equation 3)) . The scatter in Young’s approxi¬ 
mation of the scintillation noise is purely due to varying air- 
mass during the observations. The scatter in the SCIPAR 





Figure 5. Scintillation noise comparison between Youngs’s ap¬ 
proximation (equation [2]l and estimated scintillation noise from 
real SCIDAR profiles (equation IJJ. We see that Young’s approxi¬ 
mation tends to underestimate the scintillation noise as measured 
by SCIDAR. We also see that SCIDAR measures a much greater 
range of scintillation noise values. This is due to the variability 
of the atmospheric turbulence profile. The solid line represents 
ay = O'/ ij,, and is where the points should lie if equations |2] and 
|4] agree. 


scintillation noise is due to the varying airmass, the contri¬ 
bution of the wind direction component in the exponent of 
the airmass, and the inherent atmospheric variability; such 
variability is not included in Young’s approximation. 

Figure | 6 ] shows the distribution of the scintillation noise 
predicted from the SCIPAR profiles using equation |4] The 
assumed telescope size was 1 m and the exposure time was 
1 s. The normalised scintillation noise has a log-normal dis¬ 
tribution with median 0.0036 and 1st and 3rd quartiles of 
0.0028 and 0.0048 respectively. For comparison, the scintil¬ 
lation noise at zenith as estimated by Young’s equation is 
0.0026. Clearly, there is a large variability in observed scin¬ 
tillation noise which is not reflected in Young’s single-value 
approximation (equation [ 2 j . 

In figure [3 we compare measured scintillation statis- 
tics from several m ajor observatories around the world 
dKornilov et al.|[2012h including the values we have deter¬ 
mined for La Palma in figure| 6 ] We show the median, 1st and 
3rd quartile scintillation noise and the calculated Young’s es¬ 
timate for each site. From this plot we see that all of the ob¬ 
servatories in this study show similar scintillation statistics. 
However, La Palma and Manna Kea have the lowest me¬ 
dian scintillation noise. San Pedro de Martir (Mexico), Ar- 
mazones and Paranal have the highest median scintillation 
noise. Manna Kea has the smallest interquartile range sug¬ 
gesting stable low-scintillation conditions, whilst San Pedro 
de Martir has the largest interquartile range, demonstrating 
that the conditions can be more variable there. However, it 
should be noted that the data shown here is for a limited 
set of data and as we have shown scintillation noise can be 
extremely variable. In addition, we are only examining the 
variance of the intensity fluctuations and have not included 
other parameters such as the characteristic timescale, which 
will have a strong dependence on the typical wind velocity 
profile and may be important for high-speed photometry. 
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Figure 4. Atmospheric optical turbulence profiles for the night of 2014/10/09 from stereo-SCIDAR on the INT. The upper plot shows 
the turbulence strength profile and the lower plot shows the wind velocity profile. The colour-scale indicates the strength of the optical 
turbulence and the arrow length and direction indicate the turbulence speed and direction at a given time and altitude. 


Figure [7] confirms that Young’s approximation underes¬ 
timates the measured scintillation noise at all sites, typically 
by a factor of 1.5. During any one night the scintillation 
noise will usually not be equal to the median value and a 
large variation in scintillation noise can be observed, even at 
the best locations. Figure [7] also shows that the best astro¬ 
nomical sites for seeing are not necessarily the best for scin¬ 
tillation. For example. Man na Kea and Paranal have similar 
seein g statistics (for example lSarazin et ahlliooil i lChun et al.l 
l2009h but the former is better than the latter for scintillation 
noise according to the data shown in figure [7] One explana¬ 
tion for this could be that the high-altitude turbulence is 
weaker and more stable at Mauna Kea. 

It is possible to add a coefficient to Young’s approxi¬ 
mation to make it match the measured median scintillation. 
The modified Young’s approximation becomes 

(jy = 10 X 10~®C'yIl“'^^®f“^(cos7)“® exp (— 2/iobs/77), (7) 

where Cy is the empirical coefficient listed in table [T] for 
several major observatories (with a mean value of 1.5). The 
coefficients for the 1st and 3rd quartiles are also listed to 


Table 1. Values of the empirical coefficient, Cy, in our modified 
Young’s approximation (equation^ at a selection of observatories 
for the median (Cniedian)i 1st quartile (Cqi) and 3rd quartile 
(Cqs) of measured scintillation noise. 


Observatory 

^median 

Cqi 

Cq3 

Armazones 

1.61 

1.30 

2.00 

La Palma 

1.30 

1.02 

1.62 

Mauna Kea 

1.63 

1.34 

2.02 

Paranal 

1.56 

1.27 

1.90 

San Pedro de Martir 

1.67 

1.32 

2.14 

Tololo 

1.42 

1.17 

1.74 


enable the variability of the expected scintillation noise at 
each location to be estimated. 

We recommend that astronomers use our modified form 
of Young’s approximation (equation [Til with the Cy values 
listed in table [T] to give more reliable estimates of the scin¬ 
tillation noise. For even more precise estimates, or if the 
observing site is not listed in table [T] it is best to use me- 
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Figure 6. Probability Density Function (PDF) of scintillation 
noise as estimated from real SCIDAR profiles (equation for 
aim telescope and 1 s exposure times. The solid and dashed 
red lines indicate the median and 1st and 3rd quartiles, respec¬ 
tively. The solid black line indicates the value given by Young’s 
approximation (equation [2]l at the zenith. 



Figure 8. Comparison of measured noise, Cnieasuredi with esti¬ 
mated noise, ^estimated for photometric data, including scintil¬ 
lation, shot and detector noise. The measured noise values were 
obtained from light curves observed with the 4.2 m WHT (points) 
and a 0.5 m telescope on the same site (crosses), using exposure 
times between 5 and 10 s. The estimated noise values were ob¬ 
tained from SCIDAR profiles. The correlation coefficient is 0.96. 



Figure 7. Expected and measured scintillation noise for major 
observatories around the world. The bottom and top of the boxes 
indicate the 1st and 3rd quartiles and the dashed line shows the 
median measured scintillation noise. The solid line external to 
the box indicates the expected value given by Young’s approxi¬ 
mation (equation [2]l . The data for La Palma is calculated from 
recent results fro m Stereo-SCIDAR on t he INT, and the others 
are extracted from iKornilov et al.l ||2012| V 


dian atmospheric turbulence profiles (which exist for most 
observatories) in equation [I] 


3.4 Real-time scintillation estimation 

Real profiles from stereo-SCIDAR show that the optical tur¬ 
bulence profile (both turbulence strength and velocity) can 
evolve rapidly. We have seen periods of strong, high layers 
which survive for several hours and we have also seen small 
bursts of activity lastin g only a few minutes llSheoherd et al.l 
I2OI3I : lAvila et al.lll99^ . The altitude and strength of these 
layers varies with time, having significant impact on the scin¬ 
tillation noise during observations. This variability in mea¬ 
sured scintillation noise can be seen in figures |4] - [7] 

Access to concurrent turbulence profiles allows the real¬ 
time estimation of scintillation noise using equations|3]|4]and 
As the scintillation noise is dominated by high-altitude 
turbulence it is likely to be isotropic and consistent over the 
entire site. This can be validated by the fact that we see 
the same high-altitude turbulence profile before and after 
we change target with the profiler, with one target setting 
in the west as the other rises in the east. This is not true 
with low-altitude turbulence which is dependent on the local 
topography, including the telescope dome and structure. 

FigureO shows a comparison of the noise measured from 
light curves with the noise estimated from SCIDAR profiles 
obtained at the same site and time, including the effects of 
shot and detectors noise. In most cases s cintillation was seen 
to be a significant source on noise. See iFohringl ||2014|) for 
details. The large correlation coefficient of 0.96 implies that 
concurrent atmospheric optical turbulence profiles from an 
external profiler can be used to precisely model the noise 
in a photometric light curve. This enables the scintillation 
noise of an observation to be estimated, which can then be 
used for instrument characterisation, data validation and 
theoretical model fitting. 
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Figure 9. Theoretical long exposure scintillation noise as a func¬ 
tion of exposure time and telescope diameter. The scintillation 
noise was calculated for median atmospheric conditions on La 
Palma and varies between 1% for small telescopes and short ex¬ 
posure times, and 0.01% for larger telescopes and longer exposure 
times. 


4 SCINTILLATION NOISE IN PHOTOMETRIC 

MEASUREMENTS 

Figure|9]shows the expected long exposure scintillation noise 
for varying exposure times and telescope diameters. The 
scintillation noise was calculated for median atmospheric 
conditions on La Palma from equation |4] and varies between 
1 % for small telescopes and short exposure times, and 0.01% 
for larger telescopes and longer exposure times. 

The fractional shot noise {sfS/S, where S is the signal) 
can be estimated using 

as = 1/AtriS\(f), (8) 

where A is the collecting area of the telescope, rj is the 
throughput (including optics and CCD quantum efficiency), 
5\ is the bandwidth and (f> is the flux of photons incident on 
the telescope. 

When observing bright targets the sky background and 
readout noise can be considered negligible in comparison 
with the shot noise and scintillation noise. Therefore, we 
only consider the latter two noise terms. The scintillation 
noise is given by, ai = where or is the total 

photometric noise. The fractional shot noise is given by as = 
l/\/(5), where (S) is the time-averaged signal, where the 
scintillation noise is expected to approach zero. 

Figure [10] shows the estimated shot and scintillation 
noise in the short- and long-exposure regimes from on- 
sky data obtained with the 1 m Jacobus Kapteyn Tele¬ 
scope (JKT) on La Palma. As the exposure time is in¬ 
creased, the scintillation noise remains constant until we 
reach tknee ~ 0.05 s; at this point the scintillation changes 
from the short-exposure regime (with no temporal averag¬ 
ing), described by equation [3| to the long-exposure regime 
(equation |4|). For exposures longer than tknee the scintilla¬ 
tion noise reduces with a 1/v7 dependence. It is interesting 
to note that in the long-exposure regime, the fractional shot 
noise has the same time dependence as scintillation noise. 
Therefore, the ratio of scintillation noise to shot noise is 


Figure 10. Photometric noise as a function of exposure time 
on the JKT, for a magnitude 8 star observed on 21/07/2013. 
The points indicate the estimated scintillation noise, aj, and the 
crosses are the estimated fractional shot noise, Us- The lines are 
the theoretical counterparts to each segment given by equations [3] 
(short exposure; dot-dashed line), |4](long exposure; dashed line) 
and |8](fractional shot noise; dotted line). 


conserved for all exposure times. The significance of this is 
that although longer exposure times will reduce the total 
noise, the ratio of the shot noise to the scintillation noise 
will remain fixed. 

Figure ITT] shows a 2D plot of the theoretical parameter 
space for the ratio of the scintillation noise to the shot noise 
for varying telescope diameters and object magnitudes. We 
see that when observing bright objects the scintillation noise 
dominates over the shot noise, and hence dominates the total 
photometric noise. 

By combining equation [S] with equation |3] or [T] we see 
that the ratio of scintillation noise to photon noise scales 
with telescope diameter as cr/,se/o's oc and CT/,ie/crs oc 

for short and long exposures respectively. This is a 
weak dependence and so we can say that for median condi¬ 
tions and regardless of telescope diameter, scintillation will 
be greater than shot noise for (U-band) magnitudes less 
than ~13 for long exposures and ~8 for short exposures. 
It will dominate at magnitudes greater than ~12 for long 
exposures and 6.5 for short exposures. These results are for 
instruments with 100% throughput and at a zenith distance 
of 30 degrees; for other values, the results must be scaled 
accordingly. 


5 SCINTILLATION ON LARGE AND 
EXTREMELY LARGE TELESGOPES 

Kolmogorov’s model of turbulence used for the preceding 
equations assumes that energy is injected into the turbu¬ 
lent medium at large spatial scales (the outer scale, Lq) and 
forms eddies. These then break down into smaller eddies in 
a self-similar cascade until the eddies become small enough 
that the energy is dissipated by the viscous properties of the 
medium. This will occur at the inner scale, Iq. In the inertial 
range between the inner and outer scales Kolmogorov pre- 
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Figure 11. Theoretical parameter space plots for the ratio of the scintillation to shot noise in the short exposure regime (left) and the 
long exposure regime (right), for varying telescope diameter and target stellar magnitude (y-band). The short exposure time is set to 
2 ms. The long exposure time is irrelevant as both noise sources have the same exposure time dependence, making the ratio independent 
of exposure time. The black dotted line shows where the scintillation noise equals the shot noise. For any telescope diameter / target 
magnitude combinations below this line, the scintillation noise is greater than the shot noise and vice versa. The red line composed of 
circles indicates a ratio of 2, i.e. when the scintillation noise is twice the shot noise. The blue line composed of triangles indicates the 
point where the scintillation noise is an order of magnitude larger than the shot noise. 


dieted a turbulent phase power spectrum with distribution, 
^- 11/3 

As scintillation is an effect of the curvature of the wave- 
front it has power on all spatial scales up to the outer scale 
of the turbulence. On smaller telescopes the effect of the 
outer scale is negligible as we can assume the outer scale 
to be considerably larger than the telescope diameter. How¬ 
ever, as telescope sizes approach that of the outer scale (of 
the order 10 - 100 m) it should be included in theoretical 
calculations. The inner scale is generally neglected as it is 
accepted to be a few millimetres in the atmospheric case and 
too small to have any significant impact. 

All of the equations above dUElllE] and [Til assume in¬ 
finite outer scale, i.e. Kolmogorov statistics. The next gen¬ 
eration of extremely large telescopes of the scale 20 to 40 m 
are currently under construction. These telescopes will be of 
the same spatial scale, and potentially even larger, than the 
outer scale of the turbulence. If we want to examine the dy¬ 
namics of systems with spatial scales approaching the outer 
scale then we need to modify the equations that we use. 

The scintillation index for extremely large telescopes is 
expressed as (see appendix 0 for derivation) 

2 1 o 0/1 n—/ \a f Cn{h)h^ 

= 12.24L> ' t (cos 7 ) J 

j {< 1 ^ + iD/Lo{hf )) q'^{3i{Tvq)fdqdh, (9) 

where q — Df, and can be solved numerically for the par¬ 
ticular telescope and the prevailing atmospheric turbulence 
conditions. 

Figure [12] shows the expected scintillation noise for me¬ 
dian conditions on La Palma for telescope sizes between 1 m 
and 40 m in the long exposure regime. The expected scintil¬ 
lation noise is always lower than that predicted for infinite 



Figure 12. Theoretical scintillation noise for telescope diameters 
between 1 m and 40 m for the median atmospheric turbulence 
profile on La Palma in the long exposure regime (t = 1 s). The 
dashed line is the theoretical value calculated assuming an infinite 
outer scale (i.e. Kolmogorov turbulence; equation |4]l. The solid 
lines indicate finite outer scales with values of 100, 50, 20, 10, 
and 5 m. An outer scale of 10 m results in the scintillation noise 
being up to an order of magnitude less than the value given by 
equation]^ for a 39 m telescope, such as the European Extremely 
Large Telescope. 

outer scale, i.e. Kolmogorov turbulence (equation [4} , and 
this difference increases with increasing diameter. 

The factor by which equation [T] overestimates the scin¬ 
tillation noise is dependent on Lq. The outer scale of tur¬ 
bulence is currently poorly understood and estimates of its 
value are uncertain and varied. Therefore the effect of outer 
scale on the scintillation noise on large telescopes is diffi¬ 
cult to estimate, but it can be seen that, even for existing 
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large telescopes of 8-10 m diameter, the outer scale will re¬ 
duce the expected scintillation noise. The next generation 
of extremely large telescopes will all have segmented pupils 
which will complicate the situation as the gaps between the 
segments will increase the hig h frequency scintilla tion com¬ 
ponents in a non-trivial way (iDravins et al.lll998l i. In order 
to estimate this effect the aperture filter function, A{f), will 
need to be defined for each case and implemented into the 
equations above. The aperture filter function of an arbitrary 
pupil function, p{x, y) defined in a cartesian coordinate sys¬ 
tem of X and y basis, is, 

A{f)=\T(P{x,y))\P (10) 

Scintillation noise estimation on segmented or non-circular 
telescopes will therefore need to be examined on a case by 
case basis and can be done by replacing the bespoke aperture 
filter function into the derivation in appendix Scintilla¬ 
tion noise on large telescope could be overestimated by an 
order of magnitude (lKorniIovll2012ah . This suggests that ex¬ 
tremely large telescopes will be even more of a valuable facil¬ 
ity for high-precision photometry than previously thought. 


6 SCINTILLATION NOISE WITH CENTRAL 

OBSCURATIONS 

It is also important to study the impact of the central ob¬ 
scuration of the telescope on scintillation noise. Larger sec¬ 
ondary mirrors lead to more scintillation noise due to the 
smaller collecting area over which the scintillation speckles 
are spatially averaged. 

The scintillation index as a function of the telescope 
central obscuration in the short-exposure regime is given by 
(see appendix IXt 

erf = 38.44U-''/® (cos7)“®y^ Cl{h)h^ 

I (g^ + (U/Lo(h)"))"“^%" 

(Ji (vrg) - eJi (7req))^/(l - p)'^dqdh, (11) 
and in the long-exposure regime is given by, 

af,e = 12.2dD-^^h-^ (C0S7)“ / 

I (g^ + (U/Lo(h)"))'“^%" 

(Ji(7rg) - eJi(7req))^/(l - e^)'^dqdh, (12) 

where e is the ratio of the diameter of the central obscuration 
to the diameter of the primary mirror. 

Using the above equations it is possible to calculate the 
scintillation noise for any telescope diameter and in any at¬ 
mospheric conditions. However, for telescopes smaller than 
approximately 5 m the outer scale of the turbulence can be 
ignored and equations O |4] and [6] can be safely used. Most 
modern telescopes have a central obscuration of approxi¬ 
mately 30% in terms of diameter, or 10% in area. In the 
long exposure regime this makes negligible difference to the 
scintillation noise. In the short exposure regime we would 
expect the scintillation noise to be increased by a factor of 
approximately 1.2. 

In the case of a large telescope, this increase in the 


scintillation noise due to the central obscuration needs to 
be compared to the expected reduction in scintilla tion noise 
due to outer scale effects (as discussed in section^ IIKornilovI 

l2012bh . 


7 DISCUSSION 

Below we summarise the key points about scintillation noise 
that astronomers should remember when performing astro¬ 
nomical photometry: 

• Scintillation is dominated by high-altitude turbulence, 
with very little scintillation being caused by low-altitude tur¬ 
bulence. This is different to atmospheric seeing which blurs 
astronomical images. Seeing is dominated by the strongest 
layers of atmospheric turbulence, irrespective of altitude. In 
fact it is the surface layer which often dominates the seeing 
llOsborn et ahll^Oldl . Therefore, it is possible for the seeing 
to be bad whilst the scintillation noise is low and vice versa. 

• The ratio of scintillation noise to shot noise is only 

weakly dependent on telescope diameter in the short 

exposure regime and in the long-exposure regime) and 
independent of exposure time. Therefore increasing the tele¬ 
scope diameter and exposure time makes little difference to 
the ratio of shot noise to scintillation noise. 

• Young’s approximation (equation is often used by 
astronomers to estimate the scintillation noise in their light 
curves. This equation, which includes the telescope diame¬ 
ter, airmass, exposure time and altitude above sea level, is 
only an approximation to the median scintillation noise and 
is not intended to be a precise estimate. From recent Stereo- 
SCIDAR atmospheric prohles on La Palma and published 
data from iKornilov et al.l (120121 ') we have found that Young’s 
approximation consistently underestimates the fraction of 
scintillation noise by a mean value of 1.5 at a range of as¬ 
tronomical sites. 

• We have presented a modified form of Young’s approx¬ 
imation (equation [7]) that uses empirical correction coeffi¬ 
cients listed in table [T] to give more reliable estimates of the 
scintillation noise at a range of astronomical sites. 

• If even greater precision is required than provided by 
our modihed form of Young’s approximation, it is possible 
to use the median atmospheric profile for an observatory in 
conjunction with equations [3] and U) Of course the instan¬ 
taneous turbulent profile is unlikely to match the median 
atmospheric turbulence profile, and hence for the ultimate 
in precision it is necessary to use contemporaneous atmo¬ 
spheric turbulence profiles with these equations. 

• The vertical prohle of optical turbulence is very vari¬ 
able, even at the World’s premier observing locations. Fig- 
ure|4]demonstrates this for one night at the Observatorio del 
Roque de los Muchachos, La Palma. Therefore, the median 
profile can not give a precise representation of the scintil¬ 
lation noise at any particular time. Here we propose to use 
atmospheric turbulence prohling instrumentation to provide 
a better estimate of the scintillation noise in high-precision 
photometric observations. Although atmospheric turbulence 
profilers remain uncommon, with the proliferation of adap¬ 
tive optics (AO) systems they will become more widely avail¬ 
able. Indeed it is possible t o derive the require d information 
from the AO system itself (ICortes et al.|[201^ ~). 
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• A more precise knowledge of scintillation noise is im¬ 
portant to enable performance assessment, calibration and 
optimisation of photometric instrumentation. It is also use¬ 
ful when fitting models to photomet ric data (for e xample, 
extrasolar planet eclipse light curves: [Fohrind[2014l ~). and to 
help develop scintillation correction concepts. 

• As telescope diameters approach the outer scale of 
the optical turbulence, the measured scintillation noise will 
be lower than that predicted for infinite outer scale by 
Kolmogorov equations [21 El 131 El and [71 This means that 
very large and extremely large telescopes (~8 - 40 m) 
are favourable facilities for high-precision photometry. We 
present modified equations (Illl and ITTI) which should be used 
to calculate the scintillation noise with large telescopes or 
telescopes with a large secondary obscuration. 

• There are several proposed scintillation mitigation tech¬ 
niques at varying s tages of developmen t, such as Conjugate- 
Plane Photometry 10sborn_et aj JgOllf) . Tomographic wave- 
front reconstructio n (lOsborn 20151') a nd active deformable 
mirror techniques JViotto et all 1201^ '). There are also pas¬ 
sive techniques such as ‘lucky photometry’, where only data 
taken during times of low scintillation noise is used, or wave¬ 
length correction ( where one wavele ngth channel can be used 
to correct another: [Kornilo31[201l[ l. However, although the 
scintillation noise is independent of wavelength as observa¬ 
tions move away from zenith the scintillation signals become 
temporally s eparated due to the chromatic dispersion in the 
atmosphere ([Dravins et al.l[l997l l. This effect is exacerbated 
when the wind velocity is aligned with the azimuthal angle 
of the target. In addition, wavelength correction will only 
work in situations where the scintillation noise definitely 
dominates. Otherwise the noise will be made worse by the 
correction process as independent noise sources (i.e. shot 
noise) will add in quadrature. It wifi also only work in the 
large telescope {D 0.1 m) regime when the scintillation 
noise is almost independent of wavelength and it will not 
work if the variability that one wants to measure is present 
at both wavelengths. 


ACKNOWLEDGMENTS 

We are grateful to the Science and Technology Facilities 
Committee (STFC) for financial support (grant reference 
ST/J001236/1). FP7/2013-2016: The research leading to 
these results has received funding from the European Com¬ 
munity’s Seventh Framework Programme (FP7/2013-2016) 
under grant agreement number 312430 (OPTICON). The 
William Herschel Telescope, the Isaac Newton Telescope and 
the Jacobus Kapteyn Telescope are operated on the island 
of La Palma by the Isaac Newton Group in the Spanish Ob- 
servatorio del Roque de los Muchachos of the Institute de 
Astrofsica de Canarias. 


REFERENCES 

Avila R., Vernin J., Cuevas S., 1998, PASP, 110, 1106 
Brown T. M., Gilliland R. L., 1994, ARAA, 32, 37 
Chun M., Wilson R. W., Avila R., Butterley T., Aviles J., 
Wier D., Benigni S., 2009, MNRAS, 394, 1121 


Cortes A., Neichel B., Guesalaga A., Osborn J., Rigaut F., 
Guzman D., 2012, MNRAS, 427, 2089 
Dhillon V. S., Marsh T. R., Stevenson M. J., Atkinson 
D. C., Kerry P., Peacocke P. T., Vick A. J. A., Beard 
S. M., Ives D. J. Lunney D. W., McLay S. A., Tierney 

C. J., Kelly J., Littlefair S. P., Nicholson R., Pashley R., 
Harlaftis E. T., K. O., 2007, MNRAS, 378, 825 

Dravins D., Lindegren L., Mezey E., Young A. T., 1997, 
PASP, 109, 725 

Dravins D., Lindegren L., Mezey E., Young A. T., 1998, 
PASP, 110, 610 

Fohring D., 2014, PhD thesis, Durham University 
Heasley J. N., Janes K., LaBonte B., Guenther D., Mickey 

D. , Demarque P., 1996, PASP, 108, 385 

Kenyon S. L., Lawrence J. S., Ashley M. C. B., V. S. J. W., 
Tokovinin A., Fossat E., 2006, PASP, 118, 924 
Kornilov V., 2011, Monthly Notices of the Royal Astro¬ 
nomical Society, 417, 1105 

Kornilov V., 2012a, Monthly Notices of the Royal Astro¬ 
nomical Society, 425, 1549 

Kornilov V., 2012b, Monthly Notices of the Royal Astro¬ 
nomical Society, 426, 647 

Kornilov V., Sarazin M., Tokovinin A., Travouillon T., 
Voziakova O., 2012, A&A, 546, A41 
Osborn J., 2015, MNRAS, 446, 1305 

Osborn J., Wilson R. W., Butterley T., Shepherd H., 
Sarazin M., 2010, MNRAS, 406, 1405 
Osborn J., Wilson R. W., Dhillon V., Avila R., Love G. D., 
2011, MNRAS, 411, 1223 

Osborn J., Wilson R. W., Shepherd H., Butterley T., 
Dhillon V. S., Avila R., 2013, in Adaptive Optics for Ex¬ 
tremely Large Telescopes 3 Stereo scidar: Profiling atmo¬ 
spheric optical turbulence with improved altitude resolu¬ 
tion 

Roddier F., 1981, in Progress in Optics Vol. 19, The effects 
of atmospheric turbulence in optical astronomy, pp 281- 
376 

Ryan P., Sandler D., 1998, PASP, 110, 1235 
Sarazin M., Melnick J., Navarrete J., Lombardi G., 2008, 
The Messenger, 132, 11 

Shepherd H., Osborn J., Wilson R. W., Butterley T., Avila 
R., Dhillon V., Morris T. J., 2013, MNRAS, 437, 3568 
Southworth J., Hinse T. G., Joergensen U. G., Dominik M., 
Ricci D., Burgdorf M. J., Hornstrup A., Wheatley P. J., 
Anguita T., Bozza V., Novati S. C., Harpsoe K., Kjaer- 
gaard P., Liebig C., Mancini L., Masi G., Mathiasen M., 
Rahvar S., Scarpetta G., Snodgrass G., Surdej J., Thone 
G. G., Zub M., 2009, MNRAS, 396, 1023 
Tokovinin A., 2002, Applied Optics, 41, 957 
Tokovinin A., Kornilov V., 2007, MNRAS, 381, 1179 
Vernin J., Roddier F., 1973, J. Opt. Soc. Am., 63, 270 
Viotto V., Ragazzoni R., Bergomi M., Brunelli A., Dima 
M., Farinato J., Magrin D., Marafatto L., Nascimbeni V., 
Piotto G., 2012, in Proc. SPIE Mcao: Wavefront sensing 
only as a tool for high precision photometry? 

Wilson R. W., 2002, MNRAS, 337, 103 
Young A. T., 1967, AJ, 72, 747 
Young A. T., 1969, Appl. Opt., 8, 869 














Atmospheric Scintillation in Astronomical Photometry 11 


APPENDIX A: DERIVATION OF 
SCINTILLATION EQUATIONS 

We can derive the theoretical scintillatio n index as the inte- 
gral of the scintillation power spectrum llRoddierlll98ll l. 

POO 

<7? = / Wif)df (Al) 

Jo 

where W (f) is the irradiance power spectrum, given by 
dTokovinin l2002h . 

POO 

W(/) = 9.7xl0-®x4x(27r)® / cUzMf)S(zJ)A{f)fdz. 

(A2) 

S{z, f) is the Fresnel filter function to account for t he wave- 
front propagation and is given by sin^ (vrAa/^) /A^ (iRoddieil 
[1983). It is this function that gives the intensity fluctuations 
an intrinsic spatial scale of rp = '/\z. A{f) is the aperture 
filter function and is defined by 

A(/)=|P(P(*,^/))|^ (A3) 

where p(r, 6) is an arbitrary pupil function defined in carte¬ 
sian coordinate system, x and y. For a circular aperture, 
equal to one for < D^/i and ze ro elsewhere, 

A{f) = (2Ji{nDf)/{- kD/))^ llTokovininll2002l l . (j) is the fre¬ 
quency component of the refractive index power spectrum, 
for example, for Kolmogorov turbulence (f> = 

If the aperture cut-off frequency, /c oc 1/D is small 
(i.e. telescopes with D larger than a few tens of centime¬ 
tres, D r_F), such that irXzf^ -C 1, we can invoke 
the small angle approximation and S{z,f) can be written 
as (vrAa/^) /A^. Introducing the dimens ionless frequenc y 
q — Df, eauation lAll can be expressed by (lKornilovll2012al l. 

Use = 38.44D"’^^® Cl{z)z'^dz J il>{q)q^{Ji{TTq)fdq. 

(A4) 

For exposure times longer than the speckle crossing 
time, tc = DjVx, where V± is the projected windspeed 
perpendicular to the optical axis, high frequency intensity 
variations are averged out, (i.e. the scintillation temporal 
power spectrum is filtered so that only low frequency compo¬ 
nents, / < ttD/Vx affect the scintillation power). The scin¬ 
tillation weighting function with an exposure time, t, larger 
than the crossing time of the speckles is then multiplied by 
the correspondin g wind shear filter function As = D/ptV±q 
(lKornilovll2012al ). and, 

(Tie = 12.24D"'‘/® y ^vpz) j 

(A5) 

The integrals over q in the previous two equations are found 
by numerical integration to converge to 0.45 and 0.87 respec¬ 
tively. The scintillation index for short exposures becomes, 


D,se = 17-34D ^^^(cosy) ® 

POO 

/ h^Cl (h) dh, 

'0 

(A 6 ) 

and for long exposures, 



(jfje = 10 .(cos 7 )“ 

rh^cuh)df^ 

Jo Vph) ’ 

(A7) 


where h is the altitude of the turbulent layer, with h = 
2 cos ( 7 ), 7 being the zenith angle of the observation, V±{h) 
is the wind velocity profile and a is the exponent of the 


airmass. Note, that the value of the airmass exponent, a, will 
depend on the wind direction and vary between (cos 7 )”^ for 
the case when the wind is transverse to the azimuthal angle 
of the star, up to (cos 7 )”"^ in the case of a longitudinal wind 
direction. This difference comes from geometry, in the case 
where the wind direction is parallel to the azimuthal angle 
of the star the projected pupil onto a horizontal layer is 
stretched by a factor of 1 /cosy. 

For small apertures, where the aperture size is smaller 
than the spatial scale of the intensity fluctuations [D < 
rp), there is not enough spatial averaging to remove the 
dependance on the wavelength (the small angle approxi¬ 
mation on the Fresnel filter function is no longer valid, 
sin^ (ttA^/^) /A^ 7 ^ (ttXz P)^ /X^) and the scintillation in¬ 
dex can be approximated by (IDravins et al.l[l998l j. 

aj = 19.2A"''/® (cos 7 )““/® / (h) dh. (A 8 ) 

Jo 

If we replace the Kolmogorov power spectrum {(j) = 
y-ii/3^ in equationwith a Von Karmen spectrum, mod¬ 
ified to include a defined outer scale, Lq, of the form, 

</(/) = (/'+io(h)-")"“/® (A9) 

The integral in equations IA4I and IA5I now have to be eval¬ 
uated for the correct outer scale and become more compli¬ 
cated if the outer scale varies with altitude. 

The scintillation index is now be expressed as, 

afe = 12.24D-"/®t-i (cos 7 )“ J 

j {q^ + {D/Lo{hf)) g^(Ji(7rg))^dgd/i, (AlO) 

and can be solved numerically for the particular telescope 
and atmospheric parameter set in question. 

With a central obscuration t he aperture filter function 
can be expressed as (lYounelll969h . 

(All) 

where e is the ratio of the diameter of the secondary obscu¬ 
ration to the diameter of the primary mirror. 

Equation lAlII can be implemented in to the generalised 
scintillation index equation by replacing the (Ji( 7 r(;))^ in 
equation lAlOl with, 

(Ji( 7 rg) - eJi( 7 reg))^/(l - (A12) 

The final, generalised for outer scale and central obscu¬ 
ration, long exposure scintillation index is then given by, 

2 1 o 0/1 r)-4/3 ,-i ( \a f Cn{h)h^ 

Uie = 12.24D 't (cosy) y 

J {< 1 '^ + {D/Lo{hf )) 

(Ji( 7 r 5 ) - eJi( 7 reg))^/(l - e^)^dqdh, (A13) 
and the generalised short exposure scintillation equation is, 

= 38.44D-'^/® (cos-f)-^ J Cl{h)h^ 

I (gV(D/Lo(h)"))'“^%" 

(Ji (tt^) - eJi (rreg))^/(1 - e^)^dqdh. (A14) 


A{f) = - 

TT 


Ji(7rD/)^ fPiTTeOf) 


TvDf 


TveDf 

































